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We have explored the hidden symmetries of a generic four-parameter nearest-neighbor spin model, 
allowed in honeycomb lattice compounds under trigonal compression. Our method utilizes a sys¬ 
tematic algorithm to identify all dual transformations of the model that map the Hamiltonian on 
itself, changing the parameters and providing exact links between different points in its parameter 
space. We have found the complete set of points of hidden SU( 2) symmetry at which seemingly 
highly anisotropic model can be mapped back on the Heisenberg model and inherits therefore its 
properties such as the presence of gapless Goldstone modes. The procedure used to search for the 
hidden symmetries is quite general and may be extended to other bond-anisotropic spin models and 
other lattices, such as the triangular, kagome, hyper-honeycomb, or harmonic-honeycomb lattices. 

We apply our findings to the honeycomb lattice iridates Na2lrC>3 and LHIrOs, and illustrate how 
they help to identify plausible values of the model parameters that are compatible with the available 
experimental data. 
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I. INTRODUCTION 

When relativistic spin-orbit coupling dominates over 
the exchange and orbital-lattice interactions, the orbital 
moment L of an ion remains unquenched and a total an¬ 
gular momentum J = S + L is formed. This was known 
to happen in compounds of late transition metal ions such 
as of cobalt (see, e.g., Ref. [if) ; however, the “cleanest” ex¬ 
amples of spin-orbit coupled magnets emerged more re¬ 
cently: these are the iridium oxides Sr 2 Ir04 and Na 2 Ir03 
with perovskite and honeycomb lattice structures, corre¬ 
spondingly. 

By construction, magnetic ordering in these systems 
necessarily involves interactions between orbital mo¬ 
ments L , in addition to a conventional Heisenberg ex¬ 
change among the spin-part of total angular momen¬ 
tum J Q. Since the L- moment, hosted _by f 2g orbital 
in a crystal, is only an “effective” one [Tl. it need not 
be conserved during the electron hoppings, thus the L- 
moment exchange interactions are generally not SU (2) 
invariant Q. Moreover, the orbital moments have a 
“shape” and hence the L-interactions are anisotropic in 
real space, too, and thus strongly frustrated even on sim¬ 
ple cubic lattices. Altogether, this results in nontrivial 
L-Hamiltonians and orderings, including e.g. noncopla- 
nar (multi-Q) states, “hidden” Goldstone modes, etc. 
I,®- Via the spin-orbit coupling, these peculiar fea¬ 
tures of orbital physics are inherited by the “pseudospin- 
J” wavefunctions and interactions |6h14|. In essence, 
frustrated nature and quantum behavior of t 2fl -orbital 
moments Ena are transferred to that of low-energy 
pseudospins J. 

Depending on the electron configuration of ions, the 
ground state pseudospin may take different values J = 
0,1/2,1,..., and a variety of magnetic Hamiltonians with 
different symmetries and diverse behavior emerge in each 


case, because of different admixture of non-Heisenberg L- 
interactions. Perhaps the most radical departure from a 
conventional magnetism is realized in compounds with 
apparently “nonmagnetic” J = 0 ions, where a compe¬ 
tition between spin-orbit and exchange interactions re¬ 
sults in a nonmagnetic-magnetic quantum phase transi¬ 
tion [uiij. 

The case of pseudospin J = 1/2 iridates is of special 
interest. This is because Sr 2 Ir04 perovskite was found 
HUP to host cuprate-like magnetism, and honeycomb 
lattice iridates A 2 Ir03 (A =Na, Li) have been suggested 
[§] as a candidate material where the Kitaev model j2§} 
physics might be realized. Following this proposal, a sub¬ 
sequent work 0 has introduced the minimal magnetic 
Hamiltonian for iridates A 2 Ir0 3 : the Kitaev-Heisenberg 
model (KH model) - a frustrated spin model with many 
attractive properties. Most importantly, its phase dia¬ 
gram contains a finite window of a quantum spin-liquid 
phase which emanates from the pure Kitaev point of 
the model with a known exact solution (25|. To reflect 
the later experimental findings in iridates, such as the 
zigzag (Na 2 Ir03 f26j - l3 ]) and spiral (Li 2 Ir03 M) type 
magnetic orderings, the initially proposed model was 
modified by including longer-range Heisenberg [28[ [3C| 
or anisotropic 31] interactions, extending the param¬ 
eter range |32|-l34j, by considering further anisotropic 
terms in the Hamiltonian IdU - PHil . or by including spa¬ 
tial anisotropy of the model parameters (43[. An alter¬ 
native picture based on itinerant approach has been also 
suggested [3]. 

Despite the extensive efforts, no consensus concerning 
the minimal model for the honeycomb lattice iridates has 
thus far been reached. A reliable microscopic derivation 
of the exchange interactions is difficult and does not lead 
to a conclusive suggestion for the minimal Hamiltonian 
and its parameters. On the experimental side, the rich- 
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est information about the underlying spin model would 
be provided by mapping momentum-resolved spin excita¬ 
tion spectrum. However, due to the lack of large enough 
monocrystals, the inelastic neutron scattering (INS) has 
been performed on powders only f28[. Another possible 
probe - resonant inelastic x-ray scattering (RIXS) - suf¬ 
fers from a small resolution at present. While it could 
be successfully applied in case of perovskite iridates [23| , 
here the limitation comes from the much smaller energy 
scale of the excitations to be studied in detail by RIXS; 
however, the overall strength of mag netic interactions in 
Na 2 lr 03 has been quantified (45ll4aj. 


Nevertheless, the experimental data collected to date 
puts rather strong constraints on the possible models. 
First, the RIXS derived magnetic energies [H, [4(| (of 
the order of 40 meV) are much higher than the order¬ 
ing temperature (~ 15 K) suggesting strong frustration. 
Second, the magnetic scattering intensity, measured by 
RIXS at zero momentum, Q = 0, is as strong as else¬ 
where in the Brillouin zone, which implies a dominance 
of anisotropic, non-Heisenberg spin interactions. Third, 
the recent resonant x-ray scattering data pH] has revealed 
nearly ideal C 3 symmetry of the spin correlations in mo¬ 
mentum space. Moreover, inelastic neutron scattering 
data [28] have indicated that a spin gap, if present, would 
be relatively small (less than 2 meV). All these observa¬ 
tions taken together imply that the dominant pseudospin 
interactions in iridates are strongly frustrated, highly 
anisotropic in spin space, and yet highly symmetric in 
real space. By very construction, all these features are 
in fact the intrinsic properties of the KH model and its 
extended versions. 


The KH model, supplemented by other C 3 symme¬ 
try allowed terms (see below), is therefore physically 
sound and plausible. However, there is a problem of its 
large parameter space (four parameters even within the 
nearest-neighbor model) resulting in complex phase di¬ 
agrams, which makes the analysis of experimental data 
and the extraction of the model parameters a difficult 
task. In such cases, clarification of the underlying sym¬ 
metry properties of the model is often of a great help. In 
general, the spin-orbital models in Mott insulators pos¬ 
sess peculiar symmetries [l4| which are rooted in the 
bond-directional nature of orbitals. In this context, a 
special four-sublattice rotation §3] within spin space has 
proved itself as an extremely useful tool in the case of 
the original two-parameter KH model [ill, 0, l47l - l50| . It 
maps the Hamiltonian on itself but changes the Hamilto¬ 
nian parameters, connecting thereby different points in 
the parameter space. Being an exact transformation, it 
transfers the complete knowledge about some point in 
the phase diagram, including the groundstate, excitation 
spectrum, response functions etc., to its partner. Based 
solely on this self-duality of the model, the entire phase 
diagram could be sketched and the deep relations be¬ 
tween the phases understood. In addition, it also reveals 
points of hidden SU ( 2 ) symmetry, where the system is 
exactly equivalent to a Heisenberg model for the rotated 


spins. Given its usefulness, it is highly desirable to find 
and analyze similar transformations for the extended ver¬ 
sions of the KH model. 

In this paper, we introduce a systematic method to de¬ 
rive dual transformations of bond-anisotropic spin Hamil¬ 
tonians and demonstrate its results and their physical im¬ 
plications in the case of honeycomb iridates adopting the 
full nearest-neighbor model [3(1, .37, 39]. We find all the 
hidden SU (2)-symmetry points of the model, the most 
peculiar one being characterized by a “vortex”-like pat¬ 
tern with a six-site unit cell, and demonstrate how the 
characteristics of the hidden Heisenberg magnet manifest 
themselves in the anisotropic situations. By identifying 
the SU (2) points we characterize all the possible gapless 
Goldstone modes that may be encountered within the 
model. This is relevant in the context of real materials 
as the spin gap was found to be well below 2meV [28l.l29l]. 
suggesting a connection to some of the SU (2) points. Fi¬ 
nally, using a self-duality of the model, we will provide 
a link between our fits of the earlier Na 2 lr 03 data [32| 
and the recent experimental observation of the ordered 
moment direction [46[. We argue that this observation 
provides a direct access to the strength of the additional 
terms “extending” the KH model, and quantify the spin 
easy axis direction in terms of this “departure” from the 
pure KH model. This allows us to suggest plausible val¬ 
ues of the model parameters that are compatible with 
the current data. While we focus here on the case of a 
honeycomb lattice as realized in Na 2 lr 03 and more re¬ 
cently in RUCI 3 [HT] , the method is general and expected 
to produce interesting results also in the context of the 
new structural families of iridates - recently synthesized 
hyper-honeycomb [HI , [53] and harmonic-honeycomb lat¬ 
tices [54 55| , or the theoretically proposed hyperoctagon 
lattice [56f | . 

The paper is organized as follows. Sec. II introduces 
the Hamiltonian and discusses its parameters. Sections 
III and IV introduce the method, derive and discuss the 
main results of the paper - the hidden symmetries of the 
model. Sec. V and Appendix B discuss the implications 
of the results for honeycomb iridates. 


II. EXTENDED KITAEV-HEISENBERG MODEL 

We start by specifying the model Hamiltonian includ¬ 
ing all symmetry-allowed spin interactions on nearest- 
neighbor bonds. An ideal, undistorted, structure of the 
honeycomb Nah^Og plane is shown in Fig. [lfa). We will 
utilize its rotational C 3 symmetry and the three sets of 
parallel mirror planes containing the shared edges of the 
IrOg octahedra and cutting the Ir-Ir bonds into halves. 
The C 3 symmetry links the interactions for different bond 
directions while the presence of the mirror planes re¬ 
stricts the possible interactions for a given bond direc¬ 
tion. A trigonal distortion (compression or elongation 
along the Z- axis) fully preserves these symmetries so that 
our Hamiltonian applies in that case as well. Further- 
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FIG. 1. (Color online) (a) Top view of the honeycomb 
NaIr 2 C >6 plane, the definition of global X,Y,Z axes, and the 
xyz reference frame for the spin components. The X- and Y- 
directions coincide with the crystallographic a- and b-axes. 
The three bond directions of the honeycomb lattice are la¬ 
beled as a, b, and c, its two sublattices are labeled by A and 
B. (b) Two edge-shared IrC>6 octahedra of a c-bond and the 
definition of the local spin axes x, y, z (used in Eq. [TJ) . (c) Si¬ 
multaneous cyclic permutation of the Ir-Ir bond directions a, 
b, c and the spin components x, y , z when applying a C 3 
rotation to the model. 


more, recent experiments [46| indicate a nearly ideal C 3 
symmetry of the spin properties and hence suggest that 
additional terms, possibly induced by a monoclinic dis¬ 
tortion present in Na 2 lrC> 3 , can be neglected. Physically, 
this observation implies the robustness of the pseudospin 
wavefunctions against weak monoclinic distortions. 

The bond Hamiltonian is most compactly expressed in 
a local, bond-dependent, xyz reference frame for spins, 
presented in Fig. [ljb) for a c-bond. Due to the mirror 
symmetry, the in-bond S x component is forbidden to in¬ 
teract with the and S z components [36]. Following 
the notation of Ref. [36] we arrange the allowed terms into 
the form 

H {lj) \ lc = JS l -S j + KSfS z 

+ D (, SfSf - SfS f) + C (Sfsj + S*Sf) . ( 1 ) 

This four-parameter Hamiltonian extends the KH model 
(J- and K- terms) by the D-terrn bringing further 
anisotropy among the diagonal components of the inter¬ 
action, and the C-term determining the only symmetry- 
allowed non-diagonal element in the exchange interaction 
tensor. Parameter C would vanish for an isolated pair of 
undistorted octahedra; it becomes finite due to a trigonal 
distortion and/or due to the extended nature of orbitals 
in a crystal (“recognizing” the fact that the octahedra 


are canted relative to the crystal axis Z). 

To capture the C 3 symmetry, it is convenient to switch 
to cubic axes xyz , introduced in Fig. QJa) and pointing 
from an Ir ion to neighboring O ion positions in an ideal 
structure. The c-bond Hamiltonian in the cubic reference 
frame, as derived in Ref. |37|, reads then as: 

H {ij> ^ 0 = JS i -S j +KSfS] 

+r(s*s?+svs?)+r'(s*s*+s*s*+s?s]+s*s?), 

( 2 ) 

with the correspondence T = — D and T' = -bf 1 , often 
used below. For the other bond directions, the Hamilto¬ 
nian is obtained by a cyclic permutation [see Fig. [2(c)], 
resulting in one-to-one correspondence between the three 
types of bonds and interactions, as required by C 3 sym¬ 
metry. Physically, each type of bonds favors its own dis¬ 
tinct “orbital setup” to optimize the hopping energy, and 
this is fingerprinted in pseudospin interactions via spin- 
orbit coupling. For completeness, the Appendix A shows 
the Hamiltonian in the global axes XYZ ; it has certain 
advantages moving the bond dependence from the oper¬ 
ator forms to the coupling constants. 

Few comments are in order concerning the model pa¬ 
rameters. In general, calculation of exchange integrals in 
transition metal compounds with 90° d-p-d bonding ge¬ 
ometry is an intricate task, because more hopping path¬ 
ways are allowed as compared to a simpler case of 180° 
d-p-d bonding in perovskites (where theory f| has cor¬ 
rectly predicted the strength of dominant exchange con¬ 
stants). For instance, f 2 g orbitals may also overlap di¬ 
rectly, in addition to oxygen mediated hoppings; there is 
a large overlap between orbitals of t 2 g and e g symmetries 
(forbidden in perovskites), etc., resulting in a number 
of competing ferromagnetic and antiferromagnetic con¬ 
tributions which are difficult to evaluate, in particular in 
compounds with small Mott and/or charge-transfer exci¬ 
tation gaps. The uncertainties in interaction parameters 
U and Jh further affect the theoretical estimates. 

Initial consideration Jg] of the pseudospin one-half ex¬ 
change interactions in 90°-bonding geometry resulted in 
K = —2J (hitting a “hidden” SU(2) point by chance) 
in the cubic limit; later work [§j, !|, [HI using different 
approximations has changed this estimate both in terms 
of the signs and values of J and K, illustrating the diffi¬ 
culties described above. It was also found that the non¬ 
diagonal element Y allowed in cubic symmetry may take 
sizable values !36| 37], 39i, TCl- Further, T' is expected to 
become as large as the other parameters if trigonal split¬ 
ting A of the t 2g orbital level, caused by a compression 
along Z- axis, becomes comparable to spin-orbit coupling 
A; also, the trigonal field suppresses the parameter K. 
These trends are easy to understand: large trigonal field 
suppresses the in-plane components of orbital moment 
Lx and Ly, leaving the axial Ly component the only 
unquenched one; thus the pseudospin one-half Hamilto¬ 
nian, written most conveniently in global axes in this 
limit, may not contain anything but XX + YY and ZZ 
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type terms: Jxy (Sf Sf+Sf S'- ) + JzSfS , identical for 
all bonds. This is what has indeed been found by explicit 
calculations l(i. [do in the limit of A > A. This implies 
K = 0 and T' = T in this limit (see also Appendix A), 
while Jxy = (J — T) and Jz = ( J + 2T) may take any 
values depending on the microscopic details. Although 
this limit is not very realistic for Ir 4+ ion with large spin- 
orbit constant A ~ 0.4 eV ana, we may expect sizable 
values of both T and T' in Na 2 lr 03 where A seems to 
exceed 0.1 eV jKl,|5Sj. The role of T and V terms should 
further increase in other compounds based on pseudospin 
one-half Co 4+ , Ru 3+ . and Rh 4+ ions with smaller A. 

In general, the high-energy behavior of spins and or¬ 
bitals in transition metal compounds is well captured by 
the Kugel-Khomskii models [dj and their descendants [ 6 ]. 
However, the low-energy physics and ultimate magnetic 
“fixed-point” are heavily influenced by many “unpleas¬ 
ant” details originating from orbital-lattice coupling and 
distortions, unavoidable in real materials. In perovskites, 
the Kugel-Khomskii energy scale is given by At 2 /U in¬ 
dependent on spin-orbit coupling; however, this leading 
term drops out for pseudospins one-half in the edge- 
shared, 90° -bonding geometry i! , so the “high-energy” 
scale is set up by the subleading terms. In iridates, the 
hope [ 9 J is that the Kitaev-type coupling is the leading 
one among these subleading terms. Since this coupling is 
itself a correction to 4 t 2 /U, this expectation may or may 
not hold in reality. 

To summarize up to now: in real materials even with 
an ideal Co symmetry, all the four exchange parameters 
may play a significant role. This motivates us to regard 
the Hamiltonian m as an effective model with arbi¬ 
trary parameters, and look for some general symmetry 
arguments that may help to identify plausible parameter 
windows in the analysis of experimental data. 


III. SYSTEMATIC CONSTRUCTION OF DUAL 
TRANSFORMATIONS 

Having fixed the model Hamiltonian, we are ready to 
explore its dual transformations. By a dual transfor¬ 
mation we mean a prescription for site-dependent rota¬ 
tions in the spin space, S, = IZiSi, which transforms a 
spin Hamiltonian R(S) into a formally new Hamiltonian 
We are interested in self-dual transformations 
of H that map the model onto itself, preserving its all 
symmetry properties. That is, the rotated partner TL' 
( i ) has the same four terms albeit with different parame¬ 
ters J'K'D'C', and (ii) it respects the C 3 rotation rules 
encoded in Fig.[IJc), hence preserving the original distri¬ 
bution of the three types of bond-dependent interactions 
on a lattice. 

Starting with the JKDC Hamiltonian expressed as 
R(S) = J2(ij) SfHijSj where H,j are 3x3 matri¬ 
ces, we obtain TL(S) = TL'(S') = J2(ij) SORTS'' with 
H' rj = IZiHijlZj. For a self-dual transformation, the 


matrices Rh are identical to Hij , but the parameters 
JKDC are replaced by J'K'D'C , and the one-to-one 
correspondence between the bond directions and interac¬ 
tions remains intact. These two points in the parameter 
space are linked by the transformation and knowing the 
solution at one of the points, we may “rotate” it to the 
other one. 

In this section we give an algorithm to find the self-dual 
transformations for the extended KH model that map it 
onto itself. We have found a single self-dual transfor¬ 
mation JKDC -n- J'K'D'C' operating in full parameter 
space of the model; we will show it shortly below and 
return to it later when discussing experimental data. 

However, studying the hidden symmetries of the 
model, we have identified a number of restricted self-dual 
transformations that operate only in some regions of the 
parameter space, where constants J , K, D,C are all finite 
but obey certain relations, or some of them are simply 
zero. Our primary interest is in the special class of such 
transformations of the type Jo JKDC , which convert 
the Heisenberg model into the full JKDC model and 
vice versa. These transformations, to be discussed in the 
next section, reveal points of hidden SU ( 2 ) symmetry 
- by inverting the transformation the anisotropic model 
with the parameters JKDC can be exactly mapped back 
to the Heisenberg model with the exchange constant Jo- 

A. Algorithm 

A systematic search for the dual transformations seems 
to be an intricate task. Fortunately, it can be easily per¬ 
formed by computer on a finite cluster of the lattice using 
the following simple algorithm. We give it specifically for 
the case of a self-dual transformation: 

A) As a first step, we choose two rotation matrices 7 Zi, 
IZj on a selected bond ( ij). They have to preserve the 
JKDC -form given by diD, which leaves us with only a few 
choices, each having only one free angular parameter. 

B) Next, we randomly choose nonzero values of the 
initial parameters JKDC and use the relation H'- = 
IZiHijlZj together with the C 3 symmetry to determine 
the new Hamiltonian matrices for the three bond direc¬ 
tions. 

C) Knowing all the bond Hamiltonians, we may now 
determine further rotation matrices by utilizing relations 
of the type 1Z 3 = (H-j)~ 1 TZiHij and proceeding neighbor- 
by-neighbor. To fully determine the rotation matrices, 
about two thirds of the bonds need to be used. 

D) The bonds of the remaining third are used to check 
consistency, the Hamiltonian matrix determined by us¬ 
ing the rotation matrices belonging to the bond has to 
be identical to that determined in step B. If the total dif¬ 
ference on all the remaining bonds equals zero, we have 
just constructed a self-dual transformation. By scanning 
through the entire interval of the free parameter intro¬ 
duced in step A, we find all the self-dual transformations. 

The above procedure may be easily adapted to find 
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the dual transformations such as Jo -O- JKDC. In this 
case, in step A of the algorithm, we use the symmetry 
of the Heisenberg model and choose 1Z, as an identity 
matrix. The choice of the second matrix IZj is restricted 
by the requirement that H[- = 'IZ, H, = JlZj is of 

the JKDC form. 

By inspecting the rotation matrices of the cluster, we 
can identify the particular unit cell of the transformation. 
Note, that even if our cluster is smaller than this unit cell, 
we do not miss the corresponding transformation, so that 
the method is completely systematic 



J/Jo 

K/J 0 

(r=-£>)/Jo 

(r'^o/Jo 

71 

-1/3 

0 

2/3 

2/3 

71 

-1 

2 

0 

0 

71 

0 

-1 

-1 

0 

7171 

-1/9 

-2/3 

8/9 

- 4/9 

7171 

-2/3 

1 

1/3 

- 2/3 


TABLE I. Parameter values for the SU( 2) points in units of 
the exchange constant Jo of the hidden Heisenberg model. 


B. Self-duality of the extended Kitaev-Heisenberg 
model 


The systematic procedure described above has iden¬ 
tified only a single self-dual transformation JKDC o 
J' K' D'C'. This is not surprising given the complexity of 
the model. The corresponding parameter transformation 
may be written in a matrix form 


/A 

K' 


D' 

\c) 

~ 


0 

0 

Vo 
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+1 / 



(3) 


For convenience, we also give the transformation of the 
parameters JKTT' entering the Hamiltonian ([2]): 



A +1 

0 -3 
0 +| 
Vo -I 



(4) 


In terms of the spins, the transformation, labeled for fu¬ 
ture reference as 71, is simply a global 7r-rotation about 
the Z-axis defined in Fig. [TJa). The individual S x , S Y , 
and S z components transform according to 


71: (X',Y',Z') = (-X,-Y,Z) (5) 


at every site. By applying the transformation twice, we 
get an identity and the matrices in (HI are thus self¬ 
inverse. Despite its apparent triviality, this transforma¬ 
tion will play an essential role when discussing the real 
materials, see Sec. IVl below. 


IV. POINTS OF HIDDEN SU(2) SYMMETRY 

In this paragraph we find and characterize all the 
points of hidden 517(2) symmetry present in the ex¬ 
tended KH model. At these special points in the param¬ 
eter space, the anisotropic model can be mapped back 
to a Heisenberg ferromagnet or antiferromagnet. The 
SU( 2) points of the original KH model have been iden¬ 
tified ll| by virtue of the four-sublattice transformation 


introduced in Ref. 0. The corresponding ordering pat¬ 
terns on the honeycomb lattice are of stripy and zigzag 
type. A similar symmetry analysis of the KH model was 
performed for other relevant lattices [47:]. 

The extended KH model of course inherits the SU (2) 
points of the KH model and contains several new ones 
in addition. They are identified by dual transforma¬ 
tions of the type Jo O JKDC which is less general than 
JKDC -H- J'K'D'C'. Because of this, we obtain a rel¬ 
atively rich set of dual transformations characterized by 
two-, four-, and six-sublattice structure of the rotations. 
In terms of parameters, all the non-trivial SU (2) points 
are listed in Table H) We now proceed with the detailed 
description of the corresponding transformations. 


A. Summary of the SU(2) points and the 
corresponding rotations on the sublattices 

We first give a summary of the transformations as rep¬ 
resented by rotations in the real space. Each of them gen¬ 
erates an infinite number of orderings, since the ordered 
moment direction in the underlying Heisenberg model 
can be chosen arbitrarily. Figure [2] shows a few impor¬ 
tant examples. 

The simplest transformation 7/ is 7r-rotation about 
Z -axis at one of the two sublattices of the honeycomb 
lattice: 


r 2 ■■ (X', Y', Z’) = (A, Y, Z) (sublattice A) , 

(A', Y', Z') = (—A, —Y, Z) (sublattice B) . 

( 6 ) 

Its physical relevance is small due to the dominance of 
r'(= -^C) and the complete absence of K (correspond¬ 
ing to the case of strong trigonal field splitting, as ex¬ 
plained above). As a curiosity, if we choose the spins to 
lie in the honeycomb plane, 7/ converts the FM pattern 
to AF and vice versa. We may thus have an AF/FM or¬ 
dered pattern but the hidden nature revealing itself e.g. 
in the spin dynamics is that of Heisenberg FM/AF, re¬ 
spectively. 

The next transformation 71 has a four sublattice struc¬ 
ture depicted in Fig. 0(a) with 7r-rotations about cubic 
x, y, and z axes applied at sublattices 1,2, and 3, respec¬ 
tively, and no rotation involved at sublattice 4. Written 
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FIG. 2. (Color online) (a,b) Unit cells for the four- and 
six-sublattice transformations. (c,d) Stripy and zigzag pat¬ 
terns related to the FM and AF order of a hidden Heisenberg 
magnet via the four-sublattice transformation 7T The spins 
take the 2 -axis direction. (e,f) “Vortex”-like patterns gen¬ 
erated by the six-sublattice transformation T&. The spins 
are lying in the lattice plane in the case presented. The 
colors of the arrows in panels (d) and (f) indicate the sub¬ 
lattices of the hidden AF order, (g) Brillouin zones of the 
honeycomb (inner hexagon) and the completed triangular 
lattice (outer hexagon). The characteristic vectors Q a i c 
of the four-sublattice transformation and Q 1 23 of the six- 
sublattice transformation are shown in red and blue, respec¬ 
tively. (h) Bragg spots of the patterns in panels (c-f). The 
dot size is proportional to |Sq|. 


explicitly: 


: (x',y',z') = (x,—y,—z) (sublattice 1) , 

(x',y',z') = (—x,y,—z) (sublattice 2) , 


(x',y',z') = (—x,—y,z) (sublattice 3) , 
(x',y',z') = (x,y,z) (sublattice 4) . (7) 

This transformation, introduced earlier in Ref.0, is a self¬ 
dual transformation of the original two-parameter KH 
model and has been already heavily used in this context. 
Applying the transformation to an ordered Heisenberg 
FM/AF with the moments pointing along the 2 -axis, we 
get the stripy/zigzag order shown in Fig.[2](c,d). 

Perhaps the most surprising SU (2) point of the model 
is linked to the six-sublattice transformation T&. Its ro¬ 
tations are most conveniently described in the cubic co¬ 
ordinates. On the lattice sites 1, 3, and 5 [see Fig. [2Kb)] 
they correspond to cyclic permutations among the spin 
components. On the lattice sites 2, 4, and 6 the rotations 
correspond to anti-cyclic permutations which have to be 
followed by a spin-inversion. Altogether the transforma¬ 
tion can be written as 


7e : (x',y',z') = (x,y,z) 

(x',y',z') = {-y,-x 
(x',y',z r ) = (: y,z,x ) 
(■ x',y',z ') = (-x,-z 
(x',y',z') = (z, x, y) 
(x',y',z') = (- z,-y 



(sublattice 1) , 

-*) 

(sublattice 2) , 


(sublattice 3) , 

- y ) 

(sublattice 4) , 


(sublattice 5) , 

-x) 

(sublattice 6) . (8) 


It is easy to see, that for I\ = F(= —D) and J = T' = 0, 
these rotations lead to the isotropic Heisenberg Hamil¬ 
tonian. As an example, we consider the c-bond 1-2 of 
Fig. (21(b). By exchanging x and y at site 2, the non¬ 
diagonal T-term in fl2J) becomes diagonal and the inver¬ 
sion ensures its proper sign. Sample patterns generated 
by Tq and showing a “vortex”-like structure are presented 
in Fig. EKe,f). The peculiarity of the 517(2) points is 
now best demonstrated: the Hamiltonian is completely 
anisotropic containing K and T(= —D) terms only, the 
ordered spins form a very unusual pattern, yet the hidden 
nature of the system is exactly that of the Heisenberg FM 
or AF, including e.g. the presence of gapless Goldstone 
modes. 

Apart from revealing a hidden 517(2) point of the 
present model, the Tq transformation has a remarkable 
property that deserves a special attention. Namely, ap¬ 
plying 7(3 to the Kitaev Hamiltonian, we notice that it 
re-distributes three types of Ising-interactions on a hon¬ 
eycomb lattice such that at each hexagon a Kekule type 
pattern is formed [6lj. We thus arrive at the so-called 
Kekule-Kitaev model [62]. In other words, the Kitaev 
and Kekule-Kitaev models are exact dual partners linked 
via the Tq transformation. This observation should be 
helpful in studying both models, in particular of their 
extended versions including a Heisenberg term [62|, [63]. 

Two more transformations providing SU (2) points are 
obtained as the combinations T 1 T 4 and T^Te- They share 
the sublattice structure with Ta and Tq, respectively. 
Adopting the extended KH model, the former one is 
probably the 51/(2) point closest to the real situation 
in Na 2 lrC >3 as will be discussed in Sec. El 
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r>o 

A: TW0.894 T 2 
B: r’= -0.398 %% 
C: r’= -0.535 %% 


□ r’= -0.365 Na 2 lr0 3 
o r’= -0.395 Li 2 lr0 3 


r<o 

A: r’= -0.894 % 

B: r’=+0.398 %% 
C: TW0.535 T 2 T 6 


FIG. 3. (Color online) (a) Depiction of the 51/(2) points 
using the parametrization of Ref. l37l . J = sin 8 cos <j>, K = 
sin 8 sin rj>, T = cos 8. The distance from the center of the circle 
corresponds to 9 going from 0 (center) through 7r/4 (dashed 
circle) to 7t/2 (solid circle). The polar angle is <j>. Filled 
squares show the 57/(2) points with T' = 0, open squares 
those with nonzero T' values given on the right along with 
the transformation label. The color of the points indicates 
their hidden FM (blue) or AF (red) nature. The green square 
(circle) shows the parameter values specified in Sec. IVl when 
discussing Na2lr0 3 (Li2lr0 3 ). (b) The same as in panel (a) 
but with r = -cos 6. 


B. Implications for the phase diagram 


After examining the nature of the individual SU (2) 
points, we want to visualize now their positions in the 
parameter space, get a sketch of the phase diagram, and 
infer the relations between the individual phases. The 
result can be compared with the published phase dia¬ 
grams of Refs. [13 and HU, obtained by classical analysis 
and partly complemented by exact diagonalization. For 
this reason, we adopt the representation of the parame¬ 
ter space introduced in Ref. 1371 The overall energy scale 
irrelevant for the phase diagram is removed and J , K, 
T are parametrized using “spherical” angles 9 and 0 via 
J = sin 9 cos </>, K = sin 9 sin cf>, and T = —D = ±cos 9, 
keeping T'(= -^=(7) as a separate parameter of the phase 
portrait. 


Shown in Fig. [3] is the complete set of SU (2) points 
of the extended KH model. The outer rings correspond 
to the original KH model and contain the trivial SU (2) 
points and the two well-known 71 hidden SU( 2) points 
of the KH model characterized by a stripy and zigzag 
pattern. Still within the JKT plane is the “vortex” 7e> 
point associated with a “vortex”-like pattern. The corre¬ 
sponding phases determined by these SU(2) points can 
be observed in Figs. 2 and 3 of Ref. H3, with the 7! point 
lying in their 120° phase. 

Three more SU(2) points A , B 1 and C characterized 
by a nonzero value of F 7 are shown as projected onto the 
JKT plane. For the T > 0 case presented in Fig. (31(a), 
they are of AF character, one of them appears for pos¬ 
itive (point A) and two for negative (points B and C) 
values of F'. The point A (given by T 2 ) of hidden AF na¬ 
ture can possess FM pattern as discussed in the previous 
paragraph. The region between the true FM Heisenberg 
point and the point A in the phase diagram obtained 
classically is therefore filled by the FM phase extending 
as T' increases [see panels (c) and (e) of Fig. 2 of Ref. l38j . 
However, the (hidden) nature of this phase changes from 
FM to AF which should manifest itself e.g. on the char¬ 
acter of the magnon dispersion. Similarly, the presence 
of the points B (7171) and C {TzTe) of zigzag and “vor¬ 
tex” character, respectively, explains the enlarged region 
of the corresponding phases in the classical phase dia¬ 
gram for T' < 0 [see the panels (a) and (d) of Fig. 2 of 
Ref. HU. We also observe an intimate relation between 
the zigzag phase emanating from the B (7171) point and 
that connected to the zigzag SU(2) point of the origi¬ 
nal KH model (given by 71). Due to the additional 71 
rotation, their ordered moment directions are related by 
7r-rotation about the global Z- axis. This point will be 
further discussed in Sec. [V] Finally, similar conclusion as 
for the r > 0 case presented in Fig. [3(a) can be drawn 
for the r < 0 case shown in Fig. [3(b). The SU{2) points 
are related by inversion with respect to the center of the 
circle and the opposite FM/AF nature. 

In summary, we have illustrated that the gross features 
of the phase diagram of the extended, four-parameter KH 
model can be deduced solely by inspecting the nature of 
the points of hidden SU (2) symmetry and their location 
in the parameter space. 


C. Spin excitation spectra 


We proceed further by inspecting the spin excitation 
spectra at the 517(2) points associated with 71 and 71 
transformations, and see how they are related to those 
of the simple Heisenberg model. To this end, the dual 
transformations have to be expressed in Fourier space 
and relations between the Fourier components S q of the 
dual partners have to be established. The situation is 
somewhat complicated by the two-sublattice structure of 
the honeycomb lattice, requiring us to introduce an ad¬ 
ditional index [see the labels A and B in Fig. UK a) for the 












convention used below]. 

In both cases, it is convenient to use the cubic axes xyz. 
The four-sublattice transformation has three character¬ 
istic vectors Qa/b = — 7t/3) and Q c = ( 0 , 27t/3) 

touching the Brillouin zone boundary in the middle of 
its edges [see Fig. (2Tg)]- The rotation matrices have a 
simple diagonal form, reflecting only the sign changes of 
the respective components 

7 Z A/B = diag (±e l Q« R ,±e lQ >’ R ,e lQ c R ) . (9) 

The six-sublattice transformation written in Fourier rep¬ 
resentation has a full matrix structure 

n A/B =±\(l + M A/Bl + M* A/B 7 *) ( 10 ) 


with the factor 7 = | + e l< ^ 2 R + e lCi3 ' R ) and the 

matrices 


I = 

(l 1 1\ 
1 11 

m a = 

( 1 c* c N 
c 1 c* 

M b = 

( c 1 c*\ 
1 c* c 


V 1 V 


c l ) 


c V 


(n) 

where c = e 2 ™/ 3 . The characteristic vectors Qi 2 = 
(—27t/3v / 3, ±27t/3) and Q 3 = (47r/3x/3,0) shown in 


Fig. [2][g) again touch the boundary of the Brillouin zone, 
now in its corners. The dual transformation takes a gen¬ 
eral form S' AR = e l Q' R lZAQS A R (here for sublattice 
A) which translates into 

S 'm = TZ-aqSA, q-Q , (12) 

Q 

i.e., the Fourier components get shifted by the charac¬ 
teristic vectors. As a side result, the above relation gives 
the Bragg spots derived from the Bragg spots of Heisen¬ 
berg FM/AF (S A:q= 0 = ±S B , q =o = 1) and presented in 
Fig- Oh). 

To study the spin excitations, we employ the spin sus¬ 
ceptibility tensor defined as 

pOO 

Xa p(q,u)=i (lS«(t),S 0 _ q m)e i{u+lS)t dt. (13) 
Jo 

It is evaluated at the SU( 2) points by first decomposing 
S q into the A and B-sublattice contributions via 

S q = e iV3 «-/ 2 (e iq v' 2 S Aq + e~ iq y/ 2 S B q) , (14) 

applying the dual transformation in the Fourier form of 
Eq.QJ to get back to the underlying Heisenberg model, 
and using the spin susceptibility for the Heisenberg model 
obtained within linear-spin-wave (LSW) approximation. 
In the case of 71, this brings simple q-shifts by Q a , Q b , 
and Q c for the individual components. For 7e>, the cor¬ 
responding expressions are somewhat more involved con¬ 
taining a non-shifted contribution and shifted contribu¬ 
tions combining pairs of the characteristic vectors Q 3 , 



FIG. 4. (Color online) (a) LSW dispersion of the Heisen¬ 
berg FM (blue) and AF (red) on the honeycomb lattice. The 
width of the lines indicates the trace of the spin susceptibility 
tensor, Y, a XadQ, w), calculated in the LSW approximation, 
(b) The same for the stripy and zigzag state presented in 
Fig- [He, d). Energy is scaled by Jo of the hidden Heisenberg 
magnet, (c) The same for the “vortex”-like patterns presented 
in Fig.[He,f). 


Q 2 , and Q 3 . Without going into details, the presence of 
both shifted and non-shifted parts can be easily under¬ 
stood based on Eq. ITU1 

Presented in Fig. [3] are the traces of the spin suscepti¬ 
bility tensor of the Heisenberg model and the extended 
KH model at the two hidden SU{2) points under con¬ 
sideration. For completeness, we demonstrate both hid¬ 
den FM and AF case characterized by quadratic and 
linearly dispersing Goldstone modes, respectively. The 
situation is more transparent for the four-sublattice pat¬ 
terns - stripy (hidden FM) and zigzag (hidden AF) - 
since the spinwave dispersions are just shifted with the 
q=M points replacing the Goldstone points T and T' 
of the Heisenberg case. In our example, we have cho¬ 
sen z-axis as the ordered moment direction. For the 
magnons, which are in fact deviations of the ordered mo¬ 
ment in x and y directions, only Q a and Q b shifts are 
active, selecting four out of the six M-points in total. 
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The remaining two are the Bragg spots reached from T 
and r by Q c shifts active for the ordered z spin compo¬ 
nent. The Bragg spots and the Goldstone points are thus 
complementary in this case. The spin excitations asso¬ 
ciated with the six-sublattice patterns are significantly 
more complicated. They contain both shifted Goldstone 
modes [in Fig. [2(c) such a mode appears at q=K -point 
coinciding with Q :i \ and Goldstone modes at the char¬ 
acteristic momenta q = T and q = T' of the underlying 
Heisenberg model. In the latter case just the intensity of 
the modes has been transferred by the dual transforma¬ 
tion, making e.g. the linear Goldstone mode at q = T 
the most intense one in the hidden AF case. 

A similar analysis of the spin excitations as presented 
here for Ta and 7(3 SU( 2 ) points can be performed for 
the remaining SU(2) points. Due to the nature of the 
relevant transformations, no other characteristic vectors 
appear. Therefore, q = T, q = T' and their counter¬ 
parts shifted by the vectors Q a b c and Q x 2 3 entering 
the transformations 71 and 7e> constitute the entire set 
of the wavevectors of the Goldstone modes that can be 
observed within the extended KH model. 


V. APPLICATION TO THE REAL MATERIALS 

The aim of this work was to study the basic symmetry 
properties of the extended KH model - a promising spin 
Hamiltonian for the magnetism of honeycomb iridates. 
Below we illustrate how this knowledge, taken together 
with the experimental data, helps to locate the plausible 
windows in otherwise very large parameter space even for 
this nearest-neighbor (NN) model. We will show that, 
despite having only a single result, the search for full 
self-dual transformations JKDC O J'K'D'C' of the ex¬ 
tended KH model provides us with a surprisingly use¬ 
ful tool in the context of Na 2 lr 03 . This utility of 71 
emerges due to the recent observation of the magnetic 
moment direction [4f| which, as we see shortly, imposes 
an important constraint on the model parameters. This 
is because, in general, the data on magnetic easy axes 
in a crystal, along with the magnon gaps and torque 
magnetometry data, provides a direct information on the 
symmetry and strength of the anisotropy terms in spin 
Hamiltonians, and the case of Na 2 lrC >3 is of course not 
any special in this sense. 

To begin with, we recall that Na 2 lr 03 shows so-called 
zigzag order, where the spins on a and b bonds are parallel 
and form ferromagnetic chains that run along X direc¬ 
tion and couple antiferromagnetically along the K-axis. 
This relatively simple collinear magnetic structure has 
been first explained [H, [30] as due to 2 nd -NN J 2 and 
3 rd -NN J 3 Heisenberg couplings (which are often relevant 
in compounds with 90°-bonding geometry - well known 
example is quasi-one dimensional cuprates). This model 
emphasizes a geometrical frustration which is realized at 
large values of 72,3 and resolved by the C 3 symmetry 
breaking zigzag formation. 


However, as argued in the Introduction, more recent 
data [45|, [46| suggests that the origin of frustrations is 
largely related to the non-Heisenberg-type interactions 
which are bond-dependent and hence highly frustrated 
even on the level of NN-models. A minimal NN-model 
of this sort is the KH model, which has been shown {32j 
to host zigzag order in its phase diagram indeed. We fol¬ 
low this way of reasoning and explore below the extended 
version of the KH model as the basic NN-model for iri¬ 
dates. On the way, we will also see the point where the 
data may require the presence of additional terms 72,3 
too, suggesting that the both “zigzag theories” above are 
the part of a full story. 

In Ref. [32j, the available experimental data on Na 2 lr 03 
has been fitted using the two-parameter KH model, re¬ 
garding it as a phenomenological spin Hamiltonian with 
arbitrary parameters. For K = 21meV and 7 = — 4meV, 
the model was found consistent with experiments in 
terms of the type of magnetic ordering, temperature de¬ 
pendence of static magnetic susceptibility and the low 
energy spin-excitation spectrum being compared to pow¬ 
der INS. Later, RIXS experiments [4f| [46[ confirmed the 
presence of a high energy branch of spin excitations, with 
an even better agreement obtained if the LSW calculation 
of Ref. 32 is replaced by a more suitable exact diagonal- 
ization [65j |. 

However, the recent data [46j on the moment direc¬ 
tion came about as an unexpected surprise, challenging 
at first glance the above coherent description of Na 2 lrC> 3 . 
The point is that within the original two-parameter KH 
model, the zigzag order is characterized by the spins 
pointing towards one of the oxygen ions [see Fig. |T](a)] . 
This expectation is generic and guaranteed by the “order- 
from-disorder” physics [ 66 ] which typically selects one of 
the high-symmetry cubic axes as the easy one, when a 
spin Hamiltonian contains the compass-type or Kitaev- 
type bond-dependent anisotropy [O,ll3|40,|67|j indepen¬ 
dent on parameter values. The resonant magnetic x-ray 
diffraction data [|i| shows instead that the magnetic easy 
axis is in fact far away from any of the Ir -0 bond direc¬ 
tions: it is oriented “nowhere” slightly below a midpoint 
between the two, x and y , oxygen ions in Figs. [2(a) and 
0(a). This is a clear indication of the significance of the 
D and C terms in the spin Hamiltonian [ 68 ]. 

To reconcile all the data at hand using now full four- 
parameter model, we first notice that the above two easy 
axis directions - the one observed in Na 2 lr 03 and the 
one expected from the KH model as used in Ref. [32| - 
are roughly related to each other simply by a 7r-rotation 
about the Z- axis. This observation gives an immedi¬ 
ate hint how to obtain a starting parameter point when 
fitting the current data set for Na 2 lr 03 within the ex¬ 
tended, four-parameter KH model in an appealingly easy 
way, and resolve the above apparent problem with the 
moment direction. 

As discussed in Sec. IIVB1 the 74 -associated zigzag 
phase of the KH model with K > 0 is related to the 
zigzag phase connected to the SU( 2 ) point B {TiTa) of 
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(a) 




FIG. 5. (a) Pseudospin angle a relative to the XT-plane (see 
inset) as a function of the parameter r = D/(K + C/y/ 2) = 
—F /(K + F'). Dashed lines show the “magic” angle ao — 35° 
and its complement do — 55°, determined by the 2 -axis and 
xy- plane, respectively, as sketched in the inset, (b) The phase 
diagram as a function of long-range couplings J 2 = J 3 and 
anisotropy parameter r. Starting with the “bare”, 7i-derived 
values of J = 5.3 meV, I\ = —7.0 meV, D = —9.3 meV, 
and C = — 6.6 meV, we have scaled D and C simultaneously 
to vary r. To stay within the zigzag phase at the smaller 
values of r < 0.59, one needs to have finite J 2,3 couplings. 
Otherwise, the NN-only extended KH model with negative 
K < 0 switches to the incommensurate and “stripy” fll [32i . 
[69j | ground states. The inset shows the exchange bonds J 2 
and J 3 . 


Fig. Ha) via 71. Remarkably, due to the nature of 71 
- a global 7r-rotation of the magnetic moments about 
Z -axis all the aforementioned consistent results [3.2] 
of the two-parameter KH model are fully preserved if 
we apply Eq. [3] to the parameters K and J of Ref. HH 
given above, the only change is the spin easy axis be¬ 
ing rotated to the proper direction as in experiment. 
The corresponding set of parameters obtained via © is: 
J = 5.3 meV, K = -7.0 meV, D = —T = -9.3 meV, 
C = y/2T' = —6.6 meV. We would like to emphasize 
that these numbers should not be taken literally; rather, 
they fix the signs of the parameters involved and put an 
upper limit for D and C, as we explain below. 

For the representative parameters given above, the 
pseudospin makes a “magic” angle of «o — 35° from 
the XY- plane, as follows from 71 construction. This is 
slightly lower than observed mm- Now, we inspect 
how the angle a varies as a function of the anisotropy pa¬ 
rameters. The result is illustrated in Fig. H a ) and shows 
that the exact value of a heavily influences the “depar- 




C [meV] 




FIG. 6 . (a) Map of the q-dependent classical energy (per site, 
in units of meV) obtained by Luttinger-Tisza method 0 ,IH 
for the “bare”, i.e. 71-derived parameters J = 5.3 meV, 
K = —7.0 meV, D = —9.3 meV, C = — 6.6 meV rele¬ 
vant for Na 2 lr 03 . The hexagon indicates the first Brillouin 
zone, (b) The same for the parameters K = D = — 10 meV, 
J = C = 0 relevant to LHIrOs. (c) Length of the order¬ 
ing vector for varying C, keeping the other parameter val¬ 
ues unchanged. The dashed (solid) line was calculated using 
the above parameters JKD relevant to Na 2 lr 03 (LHIrOs). 
Points a and b show the C values used in panels (a) and (b), 
respectively, (d) LSW dispersions for the parameters used in 
panel (a) (left) and panel (b) (right). In the latter case, we 
have taken C ~ —1.2 meV instead of C = 0 meV to stay in 
the zigzag phase at the border to the spiral phase. 


ture” from the KH model quantified by the parameter 
r = D/(K + C/y/ 2) = —Y/{K + T'). (The correspond¬ 
ing Eqs. IB131 IB141 and IB16I for the spin direction are 
derived in the Appendix B by minimizing the classical 
energy). The “magic” angle a = «o appears at r = 0.8 - 
as obtained for the above parameters. Yet, as observed 
in Fig. H a ); by increasing a for example by 10° only, 
we already find r ~ 0.3 and get closer to the \K\ \D\ 

regime. Therefore, more detailed measurements and fits 
of the ordered spin direction are highly desirable to get 
the actual values of the parameters D and C relative 
to the Kitaev term K. Doing so, it is crucial to take 
into account the fact that the pseudospin direction and 
magnetic moment direction are not the same in general; 























































11 


while they coincide in the cubic limit, a sizable trigonal- 
field splitting might be present in Na 2 Ir 03 [53 HU- It is 
thus important to quantify this splitting by independent 
measurements. 

At this point, longer-range couplings J 2j 3 become a 
part of the full spin model for iridates, for the follo wing 
reason. As a 7i partner of the zigzag phase of Ref. [32|, 
the present NN-model with large D is well in its zigzag 
ordered state. But this is not so at smaller values of D 
( e.g. for r ~ 0.5), which are required to get the spin 
angles a ~ 40° or above, see Fig. da). Incorporating 
moderate J 2 and J 3 couplings into the model, we can 
however stabilize the zigzag phase, see Fig. (5[b), and 
hence obtain the ordered spin angles above the “magic” 
one. The values of J 2 ,3 of the order of 1 — 2 meV are in¬ 
deed suggested by ab-initio calculations [11]. This shows 
again the key importance of the experimental data on 
moment directions for quantifying the balance between 
the two zigzag-supporting mechanisms discussed above: 
based on J 2i 3 geometrical frustration, and on frustration 
driven by the non-Heisenberg nature of interactions in 
spin-orbit coupled magnets. Recent observations [46j of 
a pronounced spin-space anisotropy on one hand, and an 
“intermediate” spin direction that requires finite J 2i 3 val¬ 
ues on the other hand, suggest that both mechanisms are 
at play in Na 2 IrC> 3 . 

Altogether, the present analysis using the symmetry 
properties of the model, ta king into account the recent 
data on moment direction [461 ], as well as considering 
the role of the J 2j 3 couplings suggests a plausible win¬ 
dow in the parameter space of an effective spin model 
for Na 2 IrC> 3 : J 2i 3 < J ~ |C| < \D\ < \K\, with posi¬ 
tive (AF) Heisenberg couplings J 2j 3 and J. The leading 
anisotropy terms K < 0 and D < 0 are both negative, 
while a smaller term C may in principle take any sign. 
This parameter window is globally consistent with exper¬ 
imental observations on Na 2 Ir 03 we are aware of to date, 
and may be used as a guide in future analysis, in particu¬ 
lar once q-resolved spin response becomes available, and 
the ordered pseudospin and magnetic moment directions 
(they differ in general) are obtained and confirmed by 
independent measurements. 

Even though this general result still leaves quite a free¬ 
dom, it is of great help by fixing the signs of most rel¬ 
evant couplings and their hierarchy. This is the main 
outcome of the present theory in the context of real ma¬ 
terials. Further, we note that the Kitaev coupling I\ can 
be deduced from overall magnetic energy scale, and the 
spin and magnetic moment directions should determine 
the parameter r hence D. From a careful analysis of the 
zigzag stability condition, magnon gaps and dispersions, 
paramagnetic susceptibility data, etc., one should be able 
to quantify all the model parameters including C, J , and 

7 2j 3- 

Considering this result in the context of microscopic 
theories, we notice first that the signs of J > 0 and 
K < 0 above are consistent with the original calculations 
of these parameters for honeycomb iridates ami as well 


as with the later studies j36l-[40|. Next, we may conclude 
that a contribution from t 2g — e g hopping that favors 
pseudospin interaction with K > 0 [ 6 j is not significant 
in iridates; this is also consistent with the recent calcu¬ 
lations [3a|73. Further, the present symmetry analysis 
resolves an apparent conflict with the theoretical K < 0 
inn and the positive K > 0 that follows from the best 
data-fit using the KH model [32j: in fact, the pure KH 
model with K > 0 and the extended one with K < 0 and 
sizable D,C terms are 71-dual partners (the latter one 
being physical). 

More surprisingly, relatively large D (= — F) anisotropy 
term is required to “turn” the moment direction well 
away from the pure KH model position. A positive im¬ 
plication of this observation is that this term makes it 
much easier to stabilize the zigzag order (the pure KH 
model with large K < 0 would require large long-range 
J 2 3 couplings otherwise). In a view of the discussion in 
Sec. EH this suggests a presence of sizable trigonal field 
effects in Na 2 Ir 03 . Eventually, an unusual - out of any 
crystal symmetry axis - orientation of pseudospins [46[ 
should originate from a competition among the several 
anisotropy terms K, D, and C of different symmetry and 
physical origin. 

To conclude our discussion of Na 2 IrC> 3 : it seems that 
the extended KH model, likely further “extended” by 
moderate longer-range couplings, is indeed a good can¬ 
didate model for this compound. Even though these ex¬ 
tensions (to be still quantified by future experiments) 
reduce the chances for “pure” Kitaev-model physics in 
iridates, the model itself is highly interesting due to its 
rich internal structure and hidden symmetries that we 
have uncovered in this work. 

Motivated by the above, we further consider the case 
of Li 2 IrC> 3 . Since the data is limited here, the discus¬ 
sion will be brief and suggestive only. Due to the smaller 
Curie-Weiss temperature and more “ferromagnetic” be¬ 
havior of its spin susceptibility (74, 75j, this compound 
was located closer to SU( 2) point of the KH model 
[13. Even though the parameters K = 15.8 meV and 
J = —5.3 meV given in Ref. |32j correspond to the zigzag 
phase while Li 2 Ir 03 shows a spiral magnetic ordering 
[291 ]. these parameters can be used to get a hint of the 
direction in the parameter space to consider. We there¬ 
fore transform the above parameters using ([3]) to obtain: 
J = 1.7 meV, K = -5.3 meV, D = —T = -7.0 meV, 
C = V2T' = — 5.0meV. Representing the parameters for 
Na and Li compounds obtained via 7i transformation d3j) 
in Fig. (3] we see that both are close to the SU{ 2) point 
B (7174), with Li being closer, as expected. To approach 
the spiral state observed in Li 2 IrC> 3 , we first note that, in 
first approximation, K ~ D < 0 in both cases and that 
Li compound is characterized by a much smaller J. For 
simplicity, we set J = OrneV, assume K = D = — lOmeV 
to roughly preserve the overall energy scale, and reduce 
the parameter C associated with the trigonal distortion, 
which is expected to be much smaller in Li 2 IrC >3 with 
the bond angles being closer to 90°. The Luttinger-Tisza 
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[nun maps of the classical energy for the Na and Li case 
presented in Fig. [Da,b) confirm the zigzag and incom¬ 
mensurate magnetic ordering, respectively. For (7 = 0 
the incommensurate ordering wavevector is obtained as 
§Q C [see Fig. [U[b,c)], this would predict a magnetic 
Bragg peak in powder neutron diffraction experiments 
at a |Q |-value that could be consistent with experiments 
on powder Li 2 Ir 03 [29|. Finally, Fig. [Hid) compares the 
spin excitations obtained using LSW approximation. In 
the case of Na 2 Ir 03 , the dispersion is identical to that 
presented in Fig. 3 of Ref. (H, possessing a low and high- 
energy branches. As the parameter J is reduced, these 
two branches gradually merge, leading to a steeper dis¬ 
persion compared to Na 2 Ir 03 , which might be consistent 
with powder inelastic neutron scattering experiments on 
powder Li 2 IrC >3 The predicted dispersion is illus¬ 

trated in Fig. [ 6 jd) for a point on the boundary between 
the zigzag and the spiral phase. A further minor re¬ 
duction of C to enter the spiral phase and get the proper 
ordering vector should not affect this result dramatically, 
apart from the changes at low energies forming an “hour¬ 
glass” shape characteristic to spiral magnets (see, e.g., 
Ref.fzl). 


VI. CONCLUSIONS 

To summarize, we have analyzed non-trivial symme¬ 
tries of the extended Kitaev-Heisenberg model on the 
honeycomb lattice. As a main result, we have identified 
the complete set of points in the parameter space where 
this bond-anisotropic model can be transformed to a sim¬ 
ple Heisenberg model and is therefore characterized by 
hidden SU (2) symmetry. Such a dual transformation can 
be performed using a particular choice of sublattice rota¬ 
tions of the spins, specific for each of the SU (2) points. 
The sublattice structure of the transformations creates a 
number of ordering patterns which together with the lo¬ 
cation of the hidden SU(2) points in the parameter space 
give a good overview of the global phase diagram of the 
model. In terms of the spin excitations, the hidden SU (2) 
symmetry manifests itself by the presence of Goldstone 
modes inherited from the SU( 2) symmetric Heisenberg 
FM/AF on the honeycomb lattice. Their characteristic 
vectors and even the full spin excitation spectra are eas¬ 
ily obtained by an explicit transformation of the FM/AF 
case. 

One of the special transformations linked to the hidden 
SU(2) points reveals at the same time an exact duality 
between the Kitaev and Kekule-Kitaev models; this re¬ 
sult should be useful in theoretical studies of these and 
related models. 

We emphasize that, adopting the extended KH model, 
all the above results are necessary consequences of its 
symmetry which is in turn dictated by the underlying C 3 
symmetry of the lattice. 

Having the results of the general symmetry analysis 
at hand, we were able to find the region of the parame¬ 


ter space that is consistent with the observed properties 
of the honeycomb lattice iridates Na 2 Ir 03 and Li 2 Ir 03 . 
Further, a relation between the ordered moment direction 
and the model parameters is derived, which may help to 
quantify these parameters from future experiments. 

Finally, our method to systematically explore the hid¬ 
den symmetries is general and can be applied to other 
bond-anisotropic models as well. In the context of the 
iridate materials, the symmetry analysis of the extended 
KH model on hyper-honeycomb and harmonic honey¬ 
comb lattices is of a great interest. 
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Appendix A: XYZ form of the Hamiltonian 

The Hamiltonian expressed in terms of the spin com¬ 
ponents S x , S y , and S z , corresponding to the XYZ 
reference frame in Fig. [TJa), takes the form 

H {ij) || 7 = Jxy(S?S f + SjSj ) + JzSfsf 

+ A [c 7 (Sf Sf-Sf Sj ) - s 7 (SfSj +Sf Sf )] 

-By/2 [^(SfSf+SfSf) + s^S/Sf + S/Sj )] . 

(Al) 

Here the C 3 symmetry of the model is embodied in the 
factors c 7 = cos </> 7 and s 7 = sin </> 7 , where the angles </> 7 
are determined by the bond directions: </> 7 = 0, 4^ for 

the c, a, and b bonds, respectively. In terms of the origi¬ 
nal parameters JKTV, the exchange constants entering 
ED read as 


A = iK + §(r-r'), 

(A2) 

s = iiF-i(r-r'), 

(A3) 

Jxy = J + b — r', 

(A4) 

J z = J + A + 2T' . 

(A5) 


Note that it is the A and B terms which bring about 
the bond-directionality of the interactions, and hence 
they naturally support C 3 -symmetry breaking orderings 
such as zigzag in the present model. Physically, these 
terms arise from the exchange processes that involve the 
in-plane components of orbital momentum Lx and Ly 
which “know” the bond directions, like the orbitals do in 
the Kugel-Khomskii models. 
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It is also noticed that the A and B terms change the 
^-component of total angular momentum by ±2 and ± 1 , 
correspondingly. This is because the f 2g -orbital angular 
momentum L is not a conserved quantity in a crystal, and 
this commonly shows up in effective spin Hamiltonians 
due to the spin-orbit coupling. 

A strong trigonal field splits f 2g -level such that the low¬ 
est Kramers doublet (pseudospin) wavefunctions |f), |Z) 
become simple products of Lz = ±1 and spin | 4 _), |f) 
states, correspondingly; i.e., there will be a one-to-one 
correspondence between the real spin and pseudospin di¬ 
rections. Since the total spin is conserved during the 
hoppings, pseudospin is then conserved, too. Thus, the 
spin non-conserving terms A and B must vanish in this 
limit, which implies K —> 0 and T —> T' simultaneously. 
Physically, a strong compression along trigonal axis dic¬ 
tates that this axis becomes the “easy” (or “hard”) one 
for moments. Since this limit is not realized in iridates, 
we will not use the XYZ- form of Hamiltonian in this pa¬ 
per; however, it might be useful for pseudospin one-half 
Co 4+ , Rh 4+ , and Ru 3+ compounds where the spin-orbit 
and crystal field effects may strongly compete. 


Appendix B: Analysis of the classical energy and 
moment direction 


In this appendix we show the expressions used in the 
classical energy analysis. We first give the Hamiltonian in 
its momentum space form utilized within the Luttinger- 
Tisza method. By minimizing the classical energy in the 
zigzag phase we then find the ordered moment direction. 

Transforming the spin operators via Sar = 
J2 q e lq R SAq and similarly for the H-sublattice, we cast 
the Hamiltonian to the form 



^0 1 

e 2 \ 

( ° 

ei + e 2 

ei + l\ 


1 0 

ei +|r' 

ei + e 2 

0 

e 2 + l 


\e 2 ei 

°) 

\ er + 1 

e 2 + 1 

0 ) 


(B3) 

Here the momentum-dependent factors read as 


ei 2 = e -i|(±V 3 qx +3q y ) ^ (B4) 

Viq = l + 2cos^e-^ , (B5) 

V 2 q = cos V3q x + 2 cos : ^ 2L cos pp , (B6) 

r] 3q = e~ l3qy + 2 cos V3q x . (B7) 


In the Luttinger-Tisza method (7l|, [72), the matrices 
H q /2N s {te for q running through the Brillouin zone are 
diagonalized. The q-vector and the eigenvector corre¬ 
sponding to the minimum eigenvalue then determine the 
ordering resulting on a classical level and the minimum 
eigenvalue itself gives the classical energy per site. This 
approach relaxes the spin-length constraint which should 
be checked afterward. 

Next, we evaluate the classical energy for the zigzag 
state with the ordering vector Q = Q c = (0, 27t/ 3) and 
an arbitrary ordered moment direction given by a unit 
vector u. The corresponding zigzag pattern is captured 
by ’I'q = (+^W — \ U ) T ■ Using (IB1I) . we get for the clas¬ 
sical energy per site: 


Bciass = |(Ji - I< - 2 J 2 - 3 J 3 ) + \u t Mu (B 8 ) 
with the matrix 


M = 


( 2 K 

-r + 2 r 
V r 


-r + 2r r\ 
2 k r 

r o/ 


(B9) 


n = Y J ^\ H ^ q 

Q 


with 



(Bl) 


where the q-vectors cover the first Brillouin zone of the 
triangular lattice of R. The simplest expressions for the 
6 x 6 matrices H q of the momentum-space Hamiltonian 
m are obtained using the cubic axes x, y, z. Comple¬ 
menting the interactions in Eq. [2] by long-range ./ 2 and 
J 3 , we arrive at 


Hq — N s ite 




and 


with F q = i J 2 


r \ 0 0^ 
0 1 0 

V° 0 b 


m Q 


(B2) 


or equivalently 


M = 


( 2 K 

D + V 2 C 

V -D 


D + V2C 
2 K 
-D 



(BIO) 


The ordered moment direction can now be obtained as 
the eigenvector of M corresponding to its lowest eigen¬ 
value. However, as will be clear in a moment, it is more 
convenient to switch to the reference frame which coin¬ 
cides with the local x , y, z axes for c-bonds [see Fig.QJb)]. 
The matrix M is then transformed to 


M = 


(2K-D-V2C 

0 

V 0 


0 

2K + D + V2C 
-V 2 D 



(Bll) 


Gq — 4(^l^lq H - J^V^q) 

( l 0 o N 
0 1 0 

+ \K 

f ei 0 0\ 

0 e 2 0 


T—1 

O 

O 


O 

O 

1— 1 


which can be readily diagonalized and the angle a of the 
ordered pseudospin to the XT-plane can be found. 

As discussed in the main text, if we rotate the spins 
by 180° around the global T-axis, the observed moment 
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would come close to the z-axis. Since the latter is an 
attractive point for the two-parameter KH model 0, 
we guess that this rotation will transform the actual 
J,K,D,C Hamiltonian (K < 0, large D) for Na 2 lr 03 
into an effective K' , D', C' one, with K > 0 and small 
only D' and C' values, i.e. into a nearly two-parameter 
KH model (which guarantees that the corresponding ef¬ 
fective easy axis is close to z). We therefore first apply 
71 transformation via Eq. 0 calculate the moment direc¬ 
tion for effective J', A"', D’, C, and later make use of the 
expected smallness of the transformed D'. The first two 
steps yield an analytical expression for the angle a: 

1 ( 2V2 D' \ 

a = ao + — arctan --=— (B12) 

2 \2K’ + D’ + V2C'J 

with the first contribution being the “magic” angle «o = 
arcsin(l/\/3) — 35.3° of the z-axis to the A'F-plane and 
the second contribution supposed to be small. Now, we 
return to the original spin axes by applying the 71 trans¬ 
formation again. This does not alter the angle a but 
rotates the moment into its physical position: below a 
midpoint between two oxygen ions J46J. In terms of the 
original parameters we have 


sical level (but it will turn to either x- or y-oxygen direc¬ 
tion once the order-by-disorder mechanism is switched 
on 00). The moment moves down from this posi¬ 
tion once the model is extended by D and C terms of 
a proper sign. At the parameter set given in the main 
text, the moment takes the “magic” angle. By expanding 
the arctangent near this point, we arrive at the following 
formula for the deviation from a^: 


5a = — arctan 


4^2 1 - |r 
1 +¥\ 


7 


2 v / 2 1 - | r 

7 1 + ffr 


(B14) 


with the single parameter 
D 


K 




K + F 


(B15) 


This parameter quantifies the “departure” from the KH 
model and can be measured by resonant x-ray [46j or neu¬ 
tron diffraction experiments; as mentioned in the main 
text, care has to be taken in the fits by considering the 
crystal-field effects on pseudospin wavefunctions. 

In terms of the parameter r, Eci. IB13l can be rewritten 
as 


a = ao+- arctan |^2-\/2 

For D = C = 0, this equation gives the moment direc¬ 
tion towards a midpoint of two oxygens, as expected for 
negative K values of the Kitaev coupling 0 on a clas¬ 


4A'-5D + 2vy \ 

14AT + 2317 + 7V2C J 


tan2a = 4\/2——— . (B16) 

Note, that this and the previous equations for a hold at 
finite 72,3 Heisenberg corrections as well, since the easy 
axis is determined solely by the anisotropy terms. 
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